Isospin asymmetry and type-I superconductivity in neutron star matter 
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It has been argued by Buckley et. al. Q that nuclear matter is a type-I rather than a type-II 
superconductor. The suggested mechanism is a strong interaction between neutron and proton 
Cooper pairs, which arises from an assumed U(2) symmetry of the effective potential, which is 
supposed to originate in isospin symmetry of the underlying nuclear interactions. To test this claim, 
we perform an explicit mean-field calculation of the effective potential of the Cooper pairs in a model 
with a simple four-point pairing interaction. In the neutron star context, matter is very neutron 
rich with less than 10% protons, so there is no neutron-proton pairing. We find that under these 
conditions our model shows no interaction between proton Cooper pairs and neutron Cooper pairs 
at the mean-field level. We estimate the leading contribution beyond mean field and find that it is 
is small and attractive at weak coupling. 

PACS numbers: 26.60.+c,74.20.-z,97.60.Jd 

Recently, it has been suggested that nuclear matter in neutron stars might be a typc-I superconductor. The 
astrophysical evidence is that certain neutron stars have long precession periods, and it has been suggested that 
this means the proton superconductivity cannot be type-II Q , although that inference has been contested 0, Q . 
A theoretical argument for typc-I superconductivity in neutron star cores was then presented by Buckley et. al. in 
Ref. pj. These authors assumed that the effective potential for the Cooper pair fields has a U(2) symmetry under 
rotation of proton and neutron Cooper pairs into each other. Specifically, the effective potential for the Cooper pair 
fields was assumed to take the form 

y(|A pp |MA rm | 2 )«C/(|A pp | 2 + |A„ n | 2 ) 
= -/i c (|A pp | 2 + |A„„| 2 ) + a(|A pp | 2 + |A„„| 2 ) 2 + • • • U 

The assumption that the potential is approximately a function of |A pp | 2 + |A nTI | 2 leads to |A pp | 2 |A„ n | 2 cross terms, 
which provide a strong repulsive interaction between the neutron and proton Cooper pair condensates. It is argued 
in Ref. 0] that this leads to a long-range attraction between proton flux tubes, i.e. type-I superconductivity. 

The assumed C/(2)symmetry was justified by invoking the isospin symmetry of the underlying nuclear interaction. 
The authors of Ref. [lj admit that isospin is severely broken by the constraint of electrical neutrality, which, combined 
with beta-equilibration, greatly suppresses the proton Fermi momentum relative to the neutron Fermi momentum, 
but they claim that this does not affect the interaction between the p-p and n-n Cooper pairs. This seems implausible, 
since it is well known in the theory of superconductivity that the coefficients of the terms in the Landau-Ginzburg 
effective theory, including the quartic coupling a, depend strongly on the Fermi momenta of the underlying fcrmions 
For example, for a one component Fermi gas, just below the BCS critical point, the effective Landau-Ginzburg 
potential is 

n|A| 2 ) = - M |A| 2 + a|A| 4 + ... (2) 



where /i w Np(T c — T)/T c , and Np is the density of states at the Fermi surface, so for non-relativistic fcrmions 
Np oc mpp where m is the mass and pp is the Fermi momentum. The minimum of the potential occurs at Ao = 
y/ /i/(2a), so the quartic coefficient is a — fi/ (2Aq). The binding energy of the condensate is then T^(Ao) = — ^Np Aq, 
a well-known result in the superconductivity literature We see that /i, a, and V(A ) are all very sensitive to 
the Fermi momenta of the underlying fermions. Clearly for protons and neutrons in a neutron star, which have very 
different Fermi momenta but similar pairing condensates A, this will not give an effective potential of the form 

We now back up this argument with a concrete calculation. We work with a very simple microscopic model for 
nucleon-nucleon interactions, and analyze the pairing using the mean-field approximation. We take into account the 
requirements of electrical neutrality and equilibration under the weak interactions, which disfavor neutron-proton 
pairing. At zero temperature the potential is not accurately described by an low-order polynomial like JQ), but we 
can expand around the mean field and obtain the quartic coupling between small fluctuations 0„„, </> pp , 

^ifipp, 4>nn) = h a pp (/)^ p + anntfm + a„ p 2 p 2 „ H . (3) 

By explicitly calculating these couplings we find that they do not obey the U (2) symmetry assumed in Ref. 0. In fact, 
in the mean field approximation our result (|14|) shows no interaction at all between the proton and neutron Cooper 
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pairs: a np = 0. This appears to be a generic mean-field result, and does not depend on any specific features of the 
pairing interaction. Further, and as expected from the preceding discussion, we find that a pp <C a nn . We then discuss 
the lowest-order corrections beyond the mean-field approximation by calculating the effective interaction between 
neutron and proton Cooper pairs diagrammatically. We show that this interaction is sub-leading in the coupling and 
is negligible in weak coupling. At temperatures close to the critical temperature, where Landau-Ginzburg theory can 
be used to analyze vortex structure, we find that the interaction between the pp condensate and the nn is weak and 
repulsive, not strong and attractive as Eq. Q would imply. 

Our simple model of neutrons and protons is based on the isospin-symmetric Lagrangian 
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The nucleon field N aa has isospin index a = n,p and spin index a =|,|. Repeated indices are summed. We 
immediately generalize the interaction to allow different couplings for protons and neutrons, and Fierz-transform it 
into the pairing form. We assume that all pairing is in the rotationally invariant s-wave channel, so keeping only 
those terms we obtain 



-G pp p\p\piP] - G nn n\n\nin^ 
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(5) 



The isospin-symmetric case corresponds to G pp = G nn = G np . This interaction will lead to pairing of the fcrmions at 
their Fermi surfaces, by the usual BCS mechanism. We can calculate the thermodynamic potential of the paired state 
by a standard Hubbard-Stratonovich transformation that introduces complex bosonic Cooper-pair fields A nn ,A pp , 
and A pn (for a review, see 0). The La grangian then has three parts: the kinetic term is as before, and the others 
contain the Cooper pair fields: 
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By adjusting the phases of the n and p fields we can choose A„„ and A pp to be real, leaving A pn complex. As 
is well known [8|, in neutron star matter the mismatch between the proton and neutron Fermi surfaces completely 
suppresses n-p pairing, so A pn = even at small isospin asymmetry, but for now we keep the A pn term. In the 
mean-field approximation, we neglect any space-time variation in the A fields and set them equal to their vacuum 
expectation values. We can then integrate out the fermions, giving us the volume density of the thermodynamic 
potential or grand canonical potential = E/V — /J.N/V = —p. In the rest of this paper we will loosely refer to this 
as the "thermodynamic potential" . We find 
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where we have introduced an ultraviolet cutoff A on the three-dimensional momentum integral. To obtain this 
expression, we wrote the inverse propagator in the Nambu-Gork'ov basis (pf, rif, p\, n\, pi, n^, pj, nj), and observed 
that it block-diagonalizes into two identical 4x4 blocks, 
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Note that when A p „ = 0, M further decomposes into 2x2 blocks, for (p^ , p]), (n-f, n]) etc, which describe the s-wave 
pairing of protons and neutrons respectively. In Eq. ijSJ), e p = k 2 /2m — \x p and e„ = k 2 /2m — \i n . The determinant is 
straightforward to compute and we obtain 

DetM = (cj 2 +e 2 n )(u; 2 + e 2 p ) + A 2 nn (cj 2 +e 2 p ) 
+ A^(w 2 + 4) + |A p „| 2 (u; 2 + e p e„) 

+ (A„„A PP - A 2 „/2)(A„„A PP - A; 2 /2) (9) 

The requirement that the vacuum expectation values A pp , A nn , and A prl minimize the thermodynamic potential is 
expressed in the three gap equations, 
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When A pn = 0, as in neutron star matter |£j these equations decouple into independent gap equations for A pp and 
A nn . 

Because our model interaction is so simple, it predicts pairing gaps that rise with the density of states near the 
Fermi surface, so A rises with pf- This means that it does not produce realistic neutral nuclear matter, in which the 
interaction is isospin-symmctric, and pp p <C PFn (from electrical neutrality) but nevertheless A pp ss A nn . In the real 
world, this happens because the nuclear interaction at short distance is repulsive, so even though pp n ~3> PFp, the 
neutrons end up with a similar pairing gap to the protons. In our simple model, which is strictly valid only at low 
density, the four-Fermion coupling only encodes the attractive part of the interaction through the s-wave scattering 
length. We will choose G pp to be larger than G nn so that in the neutral system the protons and neutrons have the 
same pairing gap, as in real nuclear matter. It will turn out that our essential conclusion, that there is no U(2) 
symmetry of the effective potential for the Cooper pair fields, holds irrespective of whether or not the the couplings 
Gpp and G nn are equal. It is then reasonable to guess that G np should have a value somewhere between G nn and G pp . 
To be specific, we shall employ typical values p, n ~ 60 McV (kp n ~ 335 McV) and fi p ~ 8 McV (kp p ~ 123 MeV) 
and a momentum cut-off A = 750 MeV. For these parameters the four-fermion couplings that gives A„„ ~ 1 MeV 
and A pp - 1 MeV are G nn ~ 1 x 10" 5 MeV 2 and G pp - 2 x 10~ 5 McV 2 , respectively. 

Since A„ p = in beta-equilibrated neutral nuclear matter, we can write the mean-field thermodynamic potential 




We see that this is equal to the sum of the thermodynamic potentials for two species of Cooper pair bosons that do 
not interact with each other. There are no cross-terms between A pp and A„„. The effective potential does not take 
the form JJJ . This result did not depend on the specific form of the interaction. It simply arises from the fact that 
beta-equilibrium and electrical neutrality require pf p ^ PFn, and Cooper pairing is suppressed between species with 
very different Fermi momenta. 

We may use Eq. (|12|) to investigate the nature of the quartic terms that describe the coupling between fluctuations 
in the Cooper-pair densities. Fluctuations in the pairing field <fi nn , 4>pp are defined through the following substitutions 
in Eq. I|12|) : A n „ — > A nn + (f> nn and A pp — > A pp + 4> pp where A nn and A pp are the ground state expectation values 
that satisfy the gap equations. For small fluctuations an expansion of the thermodynamic potential about the mean 
field ground state is well motivated. Retaining only the quartic terms 

^(<?W, <t>nn) = h Q-pp^p + a lm ^„ + a np 4> 2 pp <i) 2 m H . (13) 

From the preceding discussion it is clear that in the mean field approximation there are no cross-terms. The coefficients 
a nn and a pp are non zero and depend in general on the chemical potentials, A„„ and A pp . Explicitly, by Taylor 
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FIG. 1: The lowest-order mean-field contribution to a„„, the coefficient of the (f>„„ term in the effective potential, which 
describes scattering between low-momentum fluctuations <f) nn in the nn condensate. The thick lines are Nambu-Gork'ov 
neutron propagators. Each double-square vertex is an insertion of the <j> nn operator. There is a similar diagram for a pp . 



expanding Eq. (|12|l we find 
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where E n = y e n + A nn and E p = y + A^ p . These contributions can also be calculated in a diagrammatic approach 
using the Nambu-Gorkov Greens functions 0- The diagram for a nn is shown in Fig.^ an d there is an analogous one 



for a pp . We find 
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where the finite temperature Nambu-Gorkov Green's function for the neutron superfluid, expressed in the (rif,n|) 
basis (see text after Eq. ©), is 

(lW s )^ — Ef^p) \ —Ann —IU S — €„ J 

and insertions of the fluctuating pairing field are 

';;). (in 

The Matsubara frequency is w s = (2s + VjitkT . In evaluating these diagrams we ignore any momentum transfer since 
we are interested only in the low momentum fluctuations. We have explicitly checked that this gives the same result 
as Eq. I|14|) . At zero temperature and when A//i is small we obtain the following analytic expressions 
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We see that the couplings are proportional to the corresponding Fermi momenta, a nn oc kp n and a pp tx kp p , so the 
scattering of proton condensate fluctuations is much weaker than that of neutron condensate fluctuations, indicating 
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FIG. 2: The lowest-order contributions beyond mean field to a np , the coefficient of the 4>nn4>pp term in the effective potential, 
which describes scattering between low-momentum fluctuations (f>nn and <j> pp in the nn and pp condensates. The thick lines 
are Nambu-Gork'ov neutron or proton propagators. The double square vertices are insertions of the <j) nn and 4>pp operators. 
The hatched vertex is the fundamental n^np^p interaction in the Lagrangian. 



a strong breaking of any symmetry that rotates proton Cooper pairs into neutron Cooper pairs. Further, the sign of 
a nn and a pp is negative. However, since we are Taylor expanding about the global minimum of the thermodynamic 
potential, there are lower and higher order terms (not explicitly written in Eq. Q13fl ) that ensure that the system is 
stable with respect to both small and large fluctuations. This is in contrast to the usual Landau- Ginzburg analysis 
just below T c , where a Taylor expansion about the normal state is characterized by a negative quadratic coefficient 
and a positive quartic term. 

We now discuss corrections to Eq. I|12[) beyond the mean field approximation. The leading order diagram that 
contributes to the scattering between neutron and proton fluctuations is shown in Fig. EI Th° important point is that 
this diagram involves the fundamental four-fermion neutron-proton interaction G np . Evaluating the diagram, we find 
that the leading-order beyond-mean-field contribution to the effective four-point interaction between the proton and 
neutron condensate fluctuations is 

1 f d 4 k 

atrip = Gnp J Tt[Q n (k,iw)TGn(k,iu)TQ n (k,i(jj)T3\ 

/d A k 
yjTr [Q p (k, iu)rQp{k, iLo)rQ p {k, iw)r 3 ] (20) 
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where Q n (k, iuj) and Q p (k, iui) are the Nambu-Gorkov Green's functions for the neutrons and protons respectively, and 
t 3 = diag(l, —1) is the diagonal Pauli matrix 0. The dimensionless function / is denned by the integral 

M = J_ i dxxV^+T {§2+x2)5/2 - (22) 

For S <C 1 we obtain the following analytic relation 

/(5) = -l + log[5 +0[6 2 }. (23) 

To see whether beyond-mean-field corrections can raise a np to a value comparable to a nn or a pp we follow Ref. 
in denning an asymmetry parameter 

e = (a nn a pp ~ a 2 np ) / {a nn a pp ) . (24) 

(Ref. expresses e in terms of the quartic couplings a pp , a nn , a np for the expansion around A = 0, but it should 
be equally valid to expand around the mean field, which is the minimum of the thermodynamic potential.) Ref. 
concluded that type-I superconductivity requires e < 1/20. As we have seen, in the mean- field approximation e = 1. 
From Eq. (|21|l we find that the beyond-mcan-ficld corrections to e are negligible by several orders of magnitude, e is 
always within 10 -4 of unity. 



6 



All our calculations so far have been at zero temperature, in which case the expansion of the potential around the 
mean-field ground state is only valid for small deviations from that state. This makes it impossible to discuss the 
expected structure of the vortices, in which the pairing fields vary from zero to their vacuum values. In order to be 
able to say anything about the vortices, we have to work at temperatures close to the critical temperature T c , where 
a traditional Landau-Ginzburg analysis is possible, using an expansion around the zero-mean-field state with only the 
quadratic and quartic terms, 

fi(T - T c ; <p nn , <j) pp ) = -^ c nn <p 2 nn ~ ^pp^lp + a c pp^pp + a c nn cj) A nn + o° p ^ p ^„ H . (25) 

We shall now perform such an analysis. We assume that the critical temperatures of the neutrons and protons are 
not vastly different and we perturb about the Fermi gas ground state. The diagrammatic method described earlier 
permits us to calculate the coefficients in Eq. <|25|) . Since we are perturbing about the normal state, we use the Green's 
function defined in Eq. but with A„„ = A pp = 0. The quadratic coefficient 
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where (J°(fc,iw s ) = \im& n ^o Q n (k,iu s ) and Gp(k, iu> s ) = limA p -»o Gp{k, iu> s ). An explicit calculation yields 
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and A = A 2 /(2mkT). The critical temperature T c for neutron and protons is defined by the relation [i c nn = and 
Hp — 0, respectively. As expected, for T > T c we obtain /xL, > while A*^„,M™ < for T < T c . The quartic 
coefficient is 
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and a^p is given by the same relation but with the Q p (k : iuj s ) replacing G®(k,ioj s ). The Matsubara sum is easily 
evaluated and we obtain 
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where 
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The coefficient that couples between neutron and proton pairing fields is given by 
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Thus in the vicinity of the critical temperature the asymmetry parameter e can be estimated using a^ n , a pp and a c np 
calculated above. For our choice of parameters we find that T c ~ 0.6 MeV. For T = 0.4 MeV the coefficients of the 
effective potential are \i c nn = 5520 MeV 2 , fi c pp = 2010 MeV 2 , a c nn = 2627, a c pp = 968 and a c np = -3(G np /G nn ). As 
before, the asymmetry parameter e ~ 1 since a np -C a nn and a np -C a pp . The fact that a^ p is negative implies that 
the neutron(proton) superfluid density will decrease in the inner core of the proton(neutron) vortex. 

Our conclusion is that in neutral nuclear matter, the disparity between the neutron and proton Fermi momenta 
provides a strong explicit breaking of the U(2) symmetry posited in Ref. 0. This breaking is far too strong to allow 
the proposed mechanism for type-I superconductivity to operate. Although a calculation of the vortex structure is 
beyond the scope of this work, our leading order calculation of a c np provides a means to estimate the strength of the 
effective interaction between neutron and proton condensate vortices, at temperatures close to T c . We note however 
that due to the non-perturbative nature of the interaction between nucleons we cannot exclude the possibility of a 
strong coupling between the neutron and proton superfluids. In our simple model the attractive interaction between 
neutrons and protons directly leads to a negative a c np leading to a depletion of the neutron superfluid in the core of 
the proton vortex. This result is robust as long as the effective interaction between neutrons and protons is attractive 
and contradicts the predictions in Ref. Q where the neutron superfluid density increased inside the proton vortex. 
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